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Abstract 

A large number of terrestrial planets in the classical habitable zone 
of stars of different spectral types has already been discovered and 
many are expected to be discovered in near future. However, owing 
to the lack of knowledge on the atmospheric properties, the ambient 
environment of such planets are unknown. It is known that sufficient 
amount of Extreme Ultraviolet (EUV) radiation from the star can 
drive hydrodynamic outflow of hydrogen that may drag heavier species 
from the atmosphere of the planet. If the rate of mass loss is suffi¬ 
ciently high then substantial amount of volatiles would escape causing 
the planet to become uninhabitable. Considering energy-limited hy- 
drodynamical mass loss with an escape rate that causes oxygen to 
escape alongwith hydrogen, I present an upper limit for the ratio be¬ 
tween the EUV and the bolometric luminosities of stars which con¬ 
strains the habitability of planets around them. Application of the 
limit to planet-hosting stars with known EUV luminosities implies 
that many M-type of stars should not have habitable planets around 
them. 
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1 INTRODUCTION 


Two de cades after t he first confirmed discovery o f planets outside the solar 
system f Wolszczanl . 1994 ; Mayor fe Queloz . 1995 1. we know more than 1500 
confirmed planets of different mass, size and surface temperature that are or¬ 
biting aroimd stars of different spectral types. Many of the gas giant planets 
discovered are orbiting so close to their parent stars that tidal effects of the 
star and atmospheric erosion due to strong stellar irradiation pl ay do minant 


role in determining their physical properties and evolution ( hammer et al 


2003 


2008 


Baraffe et al.l.l2004t Hubbard et all 120071: lErkaev et all 120071: iPenz. Micela fe Lam men . 


Sanz-Forcada et all 20111). On the other hand, many small and possi¬ 


bly rocky planets are recently discovered that may have surface temperature 
similar to that of the Earth. Therefore, the focus has rapidly changed into 
detecting planets that may have favorable environment to harbor life. 

Classically, a habitable zone is defined as the one t hat has f a vorab le am¬ 


bient temperature to keep water in liquid state flHuand . Il959l; lHartl . 11978 ; 


Easting, Whi tm ire & Raynoldsl , 19931; iSels is et ah. 2007 ). In the recent years 


a good number of planets a re detected within t he habitabl e zone around stars 


of various spectral t ypes flUdry et all. l2007t IVogt et all 120101; iPepe et al 


2011 ; Borucki et all 20 111 2012 : Bonhls et ah . 2013 1 


It is known that ab out 75% of the stars i n the extended solar neigh¬ 
borhood are M dwarfs fjReid. Gizis fe Hawlevl . 1200211 . Analysis of data ob¬ 
tained by using the transit method indicates the occurence rate of small 
and potentially rocky planets with radius ranging from 0.5 to 2.0 R® (R® 


(Dressing’ & Charbonneau. 20131: 

KoDDarami . 

2013 

) while that by using the 

radial velocity method is 0.41 per star ( 

Bonhls et ah. 

20131. It is believed 


that a good fraction of these rocky planets are situated within the habitable 
zone of their star. Therefore, there may exist a large number of habitable 
planets in the solar neighborhood and in our galaxy. 

However, a habitable planet should have appropriate amount of hydro¬ 
gen, nitrogen, oxygen gas and water molecules that support a habitable en¬ 
vironment. Various thermal and non-thermal mechanisms cause hydrogen 
to escape the atmosphere of a terrestrial planet. Thermal mechanisms in¬ 
clude hydrodynamic escape and Jeans escape. The loss of hydrogen from a 
planetary atmosphere is limited either at the homopause by diffusion or at 
the exobase by energy. Diffusion causes substantial hydrogen loss during the 
early evolutionary period of a terrestrial planet. In today’s Earth, the atrno- 
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sphere is collisionless above the tropopause and hence the barometric laws 
break down. Therefore, Jeans law is not applicable. Jean’s law is also not 
applicable if the atmosphere is not under hydrostatic equilibrium. This may 
occur by the absorption of stellar Extreme Ultraviolet (EUV) radiation with 
wavelengths ranging from 100 to 920 A. Strong EUV irradiation heats up 
the hydrogen rich thermosphere so significantly that the internal energy of 
the gas becomes greater than the gravitational potential energy. This leads 
to the expansion of the atmosphere and powers hydrodynamic escape of hy¬ 
drogen which can drag off heavier elements if the escape flux is high enough. 
The Earth does not receive sufficiently strong EUV irradiation from the Sun 
at present. Also the total hydrogen mixing ratio in the present stratosphere 
is as small as 10~ 5 . Therefore, in today’s Earth, hydrogen escape is limited 
by diffusion through the homopause. 

The bolometric luminosity is one of the important observable stellar prop¬ 
erties that determines the distance of the Habitable Zone (HZ) or the circum- 
stellar region in which a planet can have appropriate temperature to sustain 
water in liquid state. On the other hand, the EUV luminosity determines the 
rate at which hydrogen escapes from the atmosphere of a planet in the HZ. 
If the bolometric luminosity of the star is low, the HZ is closer to the star 
and hence a planet in the HZ is exposed to stronger EUV irradiation. If the 
EUV luminosity of such faint stars is sufficiently high then a habitable planet 
would lose substantial amount of hydrogen and heavier gases. It would lead 
to oxidation of the surface and accumulation of oxygen in the atmosphere, 
scarcity of hydrocarbon and significant reduction in surface water. All these 
would make the planet uninhabitable. If a habitable planet undergoes run¬ 
way greenhouse, water vapor would reach the stratosphere where it would 
get photolyzed and subsequently should escape the atmosphere. 

In this paper, I present an analytical expression for the upper limit of the 
ratio between the EUV and the bolometric luminosities of a planet-hosting 
star of any spectral type which serves as an essential condition to ensure the 
presence of sufficient amount of hydrogen, nitrogen, oxygen, water etc. that 
supports the habitability of a rocky planet in the HZ of the star. Since both 
Leuv an d Lh are observable quantities, it will help to determine if a planet in 
the habitable zone of a star is really habitable or not. Thus the limit should 
serve as a ready reckoner for eliminating candidate habitable planets. 
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2 CRITERIA FOR PLANETARY HABIT¬ 
ABILITY 


The distance of a habitable plane t from its parent star is given by feasting. Whitmire fc Ravnoldsl . 


19931; Kopparapu et al. , 2013a) 


d 


Ls V /2 i 

LeSeff) 


( 1 ) 


where L B and L 0 are the bolometric luminosity of the star and the Sun 
respectively, dg is the distance between the Earth and the Sun (1 AU) and 
S e ff is the normalized effective stellar flux incident on the planet. 

Let Leuv be the EUV luminosity at the surface of a star. If e is the 
efficiency at which the EUV is absorbed by the planetary atmosphere then 
the sphere-averaged and efficiency corrected heating rate in the planetary 
thermosphere due to stellar EUV irradiation is 


Leuv ( e S e ffL & \ 
L b \ 4vrd| ) ' 


( 2 ) 


The heating efficiency is less than one because part of the total incident EUV 
energy drives ionization, dissociation and other reactions without striping out 
the atoms or molecules from the atmosphere. 

In order to remain habitable, S must be low enough so that the rate 
of hydrogen escape from the planetary surface due to EUV irradiation does 
not exceed a critical value that may cause the heavier constituents including 
oxygen to be entrained with the outflow of hydrogen. If the mixing ratio of 
hydrogen is low in the lower atmosphere, hydrogen escape should be limited 
by diffusion. Therefore the critical rate of energy limited loss should be 
comparable to that of diffusion limited escape under such circumstance. 

Let S c be the value of S corresponding to the critical rate of hydrogen- 
loss. Therefore, the necessary condition to prevent hydrogen loss at critical 
rate by EUV irradiation can be written as : 


Leuv 4:7rd 2 E 

L b tS e ffL Q 


( 3 ) 


Now we need to derive S c , the sphere-averaged, efficiency-corrected and 
energy-limited critical EUV flux that causes hydrogen to escape at the critical 
rate. 
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3 ENERGY-LIMITED EUV FLUX 


In order to der i ve t he energ y-limi ted c ritical EUV fl ux S r , t he formalism given 
bvlWatson. Donahue fe Walked (1198111 is adopted. IWatson. Donahue fc Walker 
(1198il l applied this formalism to estimate the mass loss of Earth’s h ydroge n 


exosp here due to solar EUV irradiation. It was also used by L amm er et al 


(120031 ) for explaining the observed extende d atmosp here for the close-in hot 
and giant transiting planet HD 209458b ( Vidal-Madiar et all. 2003 1. The 
formalism is derived from the usual steady state equations of mass momen¬ 
tum and energy conservation of a dynamically expanding, non-viscous gas 
of constant molecular weight in which the pressure is isotropic. Unlike the 
Jean’s treatment, this formalism can be applied to a dense thermosphere 
of hydrogen with a fixed temperature at the lower boundary located suffi¬ 
ciently above the homopause such that the mixing ratios of heavier gases 
are negligible as compared to that of hydrogen. All the EUV energy is ab¬ 
sorbed in a narrow region with visible optical depth less than unity and no 
EUV energy is available bellow this region due to complete absorption. The 
rate o f mass los s est imated by using this method is within a factor of few 
(IHubbard et all 120071) of the same provide d by other models inc l uding those 
that involve deta iled numeri cal solutions (Baraffe et al.2004 Tian et all 
2005 ; Yelle . 2004 ). However, Murray-Clav. Chiang fe Murray ( 2009h argued 


that for the irradiated EUV flux greater than 10 4 erg cm -2 s _1 , mass loss 
ceases to be energy limited and becomes radiation/recombin ation limit ed. 
As a consequence the formalism prescribed by Watson . Donahue & Walke r 
(119811) is not applicable in that case. 

The two equations that provide the maximum rate of hydrogen escape 
and the expansion of the atmosphere due to EUV heating are given as 
( Watson. Donahue fe Walker! . 198lh : 


£ = 


<7 + 1 


( Ai /2)( 1+g )/ 2 + 1 

Aq — Ai 


( 4 ) 


and 


P = &\ 


An — 



( 5 ) 


In 


the above equations the dimensionless parameters £, /?, Aq and Ai are 
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related with the physical parameters by 




k 2 T n 


\n§GMpTTlp[ I 


( 6 ) 



(7) 

G MpTYlH 

Ao,l T m ) 

kT 0 r 0 i 

(8) 


where r 0 = Rp is the radius of the planet where the visible optical depth 
is one and the temperature is T 0 , r\ is the radius of the region where all 
the EUV energy is absorbed. The optical depth of the atmosphere to EUV 
energy at the lower boundary tq is much greater than unity. F m is the sphere 
averaged flux of escaping particles (escape flux of particles per steradian per 
second), k is the Boltzmann constant, Mp is the total mass of the planet, 
m,H is the mass of hydrogen atom and k is the thermal conductivity of the 
gas which is parameterized by k = T q where r is the visible optical depth 
of the thermosphere such that r(A q) = 1. 


4 ANALYTICAL AND NUMERICAL SO¬ 
LUTIONS 


Equation (|3J) and equation ([5]) need to be solved numerically for an arbitrary 
value of q. However, for q — 1 we obtain analytical solutions for Ai and (3 
and hence S. Thus for q — 1, equation ((4j) gives 


2(£ 1/2 Aq - 1) 

1 1 + 2£V 2 • 

Substituting equation (J9]) in equation (J5]) , equation (JTJ) gives, 


S(q = 1) = 4£ 


V i + 2£ 1/2 ) 



( K 0 kT 2 \ 

yGMpTflp J 


(9) 


( 10 ) 


Now, S = S c when F m = F c in the expression for £ given in equation ([6]) 
where F c is the critical rate of hydrogen escape. Hydrogen diffuses before 


6 









it reaches the thermobase if the background gas of heavier species is static. 
This happens if the heavier species are not absorbed or cannot escape at the 
surface. The diffusion limit for hydrogen is achieved when the heavier gases 
attain the maximum upward velocity such that the background becomes non¬ 
static f Hunten . 19731) . Therefore if F c is greater than or equal to the diffusion 
limit, heavier specie s wo uld e scape . Th e rate o f hy drogen loss limited by 


diffusion is given by Huntcn ( 1973 ); Zahnle. Pollack fc Kasting ( 1990|) 


F H (dif fusion) = 


bg(m s - m H ) 
kT 0 (l + X s /X H ) 


( 11 ) 


where b = 4.8 x 10 17 (T / K) 0 - 75 cm 1 s 1 


is th e binary diffusion coefficient 


for the two species (IZahnle fc Kasting j , 119861) . Xp and X s are the molar 
mixing ratio at the exobase for hydrogen and the heavier atom with mass m s 
respectively and g is the surface gravity of the planet. Therefore we set 


bg(m s - m H ) 


-i-i 


F - = F ' = Rp kTj\ + A',/X H ) paltideS 51 ' S 


( 12 ) 


which gives the cross over mass m s = m H + (kT 0 F c /R 2 P ) / (bgX H ) such that 
any species with mass less than or eq ual t o m x w ould be efficiently dragged 


by the escaping hydrogen (iHunten et a.l. 119871 ). In the present derivation 


I consider m s = mo = 16m# such that the critical rate of hydrogen loss 
would enable oxygen and other atoms lighter than oxygen to be dragged 
with the escaping gas. I assume that both hydrogen and oxygen are atomic 
near t he exobase due to fast di ssoc iation of the molecules by photolysis 
( Murrav-Clav. Chiang fc Murray . 2009 ). 

Therefore, from equation (J3J) we obtain, for q = 1 


p EUV 


< 


lQnd 2 E KokTf 


C /2 (Ao« 1/2 - l) 3 


J B 


iGMpm„( 1 - e 2 )'/ 2 S e// L 0 (2^i/ 2 + 1)‘ 


(13) 


where Ao, the parameter for the planetary radius is given by equation 
and 


, _ D2 bg(m 0 - m H ) 
s — F P 


k 


1 + Xq /X E K 0 GMpmn 


(14) 


The term (1 — e 2 ) 1 / 2 in the denominator is introduced in order to include 
elliptical orbit with eccentricity e. 
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The above analytical exp ression i s der i ved b y tak ing q — 1. However, 
for a neutral gas, q ~ 0.7 ( Banks fe Kockarts . 1973f h Therefore, I solve 
equation (j4j) and equation ([5]) numerically for different values of q < 1.0 
keeping all other parameters fixed. Figure 1 presents the values of S c for 
different values of q and shows that S c (q = 1) < S c (q < 1). Therefore, 
equation (J3J) implies that the inequality given by equation ([T3]) is valid for 
any value of q < 1. It is worth mentioning that the bolometric luminosity of a 
main sequence star increases with time. For a solar-type star, the bolometric 
luminosity increases by about 10 % in every 1 Gyr. The Sun was about 30% 
fainter in the visible light during the first billion years after its birth. On the 
other hand the EUV luminosity is governed by the coronal activities of a star. 
Since a star rotates faster during younger age resulting into greater coronal 
activities, the EUV emission rate or the EUV luminosity was higher in the 
past and it decreases with time. Consequently, the rate of mass loss from 
the planetary surface was higher in the past. Since, the value of Leuv/Lb 
in the past was greater than its present value, the condition derived here is 
valid for a constant rate of mass loss corresponding to the present value of 
L E uv/Lb- A time-dependent solution would have provided more stringent 
condition. 


APPLICATION TO HABITABLE EXO¬ 
PLANETS 


The thermal c onduct ivity of h ydrog en kq = 4.45 x 10 4 ergs cm 4 s 4 K 1 


(IHanlev. McCarty fe Intermanl . 197(2) - The size and the mass of the planets 
are usually determined from the observed parameters of the planet. Transit 
method provides the size and the orbital inclination angle of the planet and 
radial velocity method provides the projected mass of the planet for known 
stellar mass. Therefore, the two unknown parameters in equation (fT3]i are 
T 0 and S eff . 

The lower boundary r 0 is located some distance above the homopause and 
the temperature T 0 at the l ower bo un dary may be fixed at the equlibrium 
temperature of the planet ( Watson. Donahue fe Walker! . 198lh . However, 
as shown in Figure 2, the above limit is weakly dependent on To because 
the parameter Ao is inversely proportional to To. Therefore, without loss 
of generality, T 0 can be fixed at 254 K, the equilibrium temperature of the 


















Earth. This makes the above expression independent of planetary albedo 
because derivation of equilibrium temperature requires the value of planetary 
albedo. With T 0 = 254K and taking Xo = 1/3 and Xh = 2/3, 1 obtain 
F c = 1.43 x 1 0 13 c m -2 s -1 w hich is two ord ers of magnitude higher than that 
calculated by Watson. Donahue fc Walker ( 198lh for the Earth irradiated by 
solar EUV radiation. This ensures the cross over mass m s is equal to mo- The 
total mass of hydrogen (Mu) i n the atmosphere. Ocean and in the crust of the 
Earth is about 1.9 x 10 23 gm flAnders fc Owen] . Il977t ISharp. Draper &; Ageel . 
20091 1. Therefore at this critical rate, an Earth like planet would lose all 


of its surface hydrogen in about 50 Myr. Water molecules would be photo- 
dissociated into hydrogen and oxygen atoms that would subsequently escape 
the planet. Even if the amount of water is not reduced sufficiently, such 
a significant loss of hydrogen and other heavier gases such as nitrogen and 
oxygen would make the environment drastically different than that of the 
pre sent Earth . 


Kopparapu et ah fl2013ai ) have presented relationships between the stel¬ 


lar effective temperature T e ff and the normalized stellar flux S e ff impinging 
on the top of the atmosphere of an Earth-like planet at the habitable zone. 
These relationships are applicable in the range 2600 K< T e ff < 7200 K. 
According to these relationships, S e ff is minimum at T e ff = 2600 K and 
increases with the increase in T e ff. S e ff also varies for different atmospheric 
conditions that determines the inner and the outer boundaries of the habit¬ 
able zone. While condition like ’’Recent Venus”, ’’Runway Greenhouse” or 
’’Moist Greenhouse” determines the inner boundary, condition such as ’’Max¬ 
imum Greenhouse” or ’’Early Mars” determines the outer boundary of the 
habitable zone. Since, Leuv /Lb hr equation (jT3j) is inversely proportional to 
S e ff, 1 use the minimum value of S e ff corresponding to T e ff = 2600K. This 
makes the upper limit on Lew/Lb independent of the effective temperature 
of the star. As S e ff increases with the increase in T eff, the value of S e ff 
corresponding to the actual T eff of the star would only tighten the limit. 
It i s wor th mentioning here that according to the relationships provided by 


Kop parapu et ah fl2013al) . S e ff increases only by 2-3 times with the increase 


in T e ff from 2600 K to 7200 K. 

Taking To = 254 K- the equilibrium temperature of the Earth, Rp = 
6.378 x 10 8 cm - ra dius o f the Earth, Mp = 5.97 2 6 x 1 0 27 gm - mass of the 
Earth, and e = 0.15 ( Watson. Donahue fc Walker! . 1981 ). I derive, from equa¬ 
tion (Jl3]) , the upper limit on L E uv/Lb for planets within the outer boundary 
of the habitable zone. The values are presented in Table 1. For the same 
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set of parameters the numerical solution of equation (J3J) and equation (ED for 
q = 0.7 is also presented in Table 1. 


Table 1: Upper limit on L EUV /L B for S e ff corresponding to T e ff = 2600K 
under atmospheric conditions that determines the outer boundary of the 
habitable zone. For q = 0.7 the results are obtained numerically. 


Conditions 

Seff 

L E ub/L b 

q=1.0 

q=0.7 

Maximum Greenhouse 

0.215 

1.33 x 10 -4 

1.56 x 10 -4 

Early Mars 

0.199 

1.44 x 10 -4 

1.69 x 10 -4 


The critical EUV flux S c is less than 10 erg cm 2 s 1 which ensures energy- 
limited escape. Also, equation ((9]) gives the expanded radius r\ = 1.57 Rp 
for T 0 = 254K. 


6 DISCUSSIONS AND CONCLUSIONS 


Estimating the stellar EUV luminosities is difficult not only because the 
energy gets completely absorbed at the uppermost layer of the Earth’s at¬ 
mosphere but also due to photoelectric absorption by the interstellar medium 
along the line of sight. Using the ROSAT space telescope, IHodgkin fe Pve 
(119941 ) estimated the EUV luminosities of a large number of nearby stars 
with spectral type ranging from F to M . Here, I use the bolometric and 
EUV lumin osities o f planet-hostin g sta rs of spectral type ranging from A to 
M, given by iSanz-Forcada et ah (20111) who derived the EUV and bolomet¬ 
ric luminosities by using the data from ROSAT, XMM-Newton and Chandra 
space telescopes. 1 have considered only those stars that are older than 1 
Gyr. It is found that out of all such stars, only two of the planet-hosting 
M stars satisfies the criteria presented in this work. The value of L EUV /L B 
for 2MASS 1207 (spectral type M8) is 1.23 x 10” 3 . For GJ 317 (M3.5), GJ 
674 (M2.5) and for GJ 176 (M2.5V), the values of L EUV /L B are 1.78 x 10 -3 , 
2.95 x 10” 4 , 4.26 x 10” 4 respectively. Therefore, none of these M stars can have 
a planet with habitable environment. GJ 832 (M1.5V) has L EB v/L B = 10“ 4 
and so it marginally satisfies the habitability criteria. However, since the 
effective temperature of this M1.5V spectral type star is higher than 2600K, 
this planet should also be considered uninhabitable. On the other hand the 
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values of L EUV /L B for GJ 436 (M2.5) and GJ 876 (M4V) are 1.58 x 10~ 5 
and 1.9 x 10~ 5 respectively. Therefore, these two M stars satisfies the EUV 
habitability criteria provided in Table 1. This is expected because M stars 
are the faintest among stars of all spectral types. As a consequence the HZ 
of M stars is located very near to the stars and hence a planet in the habit¬ 
able zone is exposed to strong EUV radiation. It is worth mentioning here 
that a rocky planet in the habitable zone of M stars may lose much of its 
volatiles during the formation because the pre-main sequence phase of such 
stars are comparatively longer. However, if the planet accretes sufficiently 
large amount of water during formation or if it were formed far a way from 
the star and then migrated to the HZ, it may remain habitable ( Lissauerl 


2007|). But strong EUV irradiation should make it uninhabitable. 


However, planet hosting stars of all other spectral types listed by lSanz-Forcada et al 


12 0111) have Leuv/Lb much lower than the upper limits provided here. The 
values of L E uv/L B for GJ 86 (K1V) is 9.77 x 10 -5 . Therefore, according 
to the present upper limit, this star marginally satisfies the habitability cri¬ 
teria of having sufficient amount of water. On the other hand, the values 
of L euv /L b for HD 87883 (KOV), e Eridani (K2V), HD 46375 (K1V), HD 
93083 (K3V), HD 130322 (KOV), HD 189733 (K1-K2) and HD 218566 (K3V) 
are 2.75 x 10“ 5 , 2.19 x 1(T 5 , 1.02 x 1(T 5 , 1.479 x 10" 5 , 1.95 x 10" 5 , 2.4 x 10" 5 
and 1.66 x 10 -5 respectively. So, all of them satisfies the upper limit for 
Leuv/L b . 

Hence, the upper limit presented here is important in deciding the hab¬ 
itability of planets around M stars. Note that the numerical values of the 
parameters involved in the present habitability condition do not differ much 
from that of the Earth. The density of a rocky planet should be about 4.0- 
5.5 gmcnr 3 . Therefore, the radius of a rocky planet does not exceed 1.6 
times the Earth’s radius and the mass must be less than 5 times the Earth’s 


mass ([Rogers I. 20151 ). More ma ssive p lanets are Neptune and Jupiter types 


of gaseous planets flMarcv et all 2014). Therefore, a large value of L EB y/L B 
for M stars may rule out the presence of habitable planets around them. 

Owing to our poor knowledge on the actual environment of exoplanets, 
the habitability of rocky planets around the habitable zone remains ambigu¬ 
ous and incomplete. However, we must impose as many conditions as can 
be determined directly and easily in order to narrow down our search for 
planets that may possibly harbor life. The present result provides such an 
important condition as it depends on two observable quantities of the planet 
hosting star - the bolometric and EUV luminosities. Spectra of habitable 
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planets around M-type stars that do not satisfy the EUV criteria presented 
here should show lack of hydrogen, nitrogen and oxygen in their atmosphere 
and hence can confirm the limit presented here. 

1 thank the reviewer for a critical reading of the manuscript and for pro¬ 
viding useful suggestions and comments. 
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Figure 1: The sphere-averaged, efficiency-corrected and energy-limited crit¬ 
ical EUV flux S c as a function of the conductivity index q. Terrestrial pa¬ 
rameters are used to calculate S c . For neutral gas, q ~ 0.7. 
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Figure 2: L EUV /L B for different values of T 0 . Solid line represents L EUV /L B 
for q = 1.0 and dashed line represents that for q = 0.7. 
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